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For Bose-Einstein condensates in double wells, iV-particle Rabi-like oscillations often seem to be 
damped. Far from being a decoherence effect, the apparent damping can indicate the emergence 
of quantum superpositions in the many-particle quantum dynamics. However, in an experiment it 
would be difficult to distinguish the apparent damping from decoherence effects. The present paper 
suggests using controlled periodic shaking to quasi-instantaneously switch the sign of an effective 
Hamiltonian, thus implementing an "echo" technique which distinguishes quantum superpositions 
from statistical mixtures. The scheme for the effective time-reversal is tested by numerically solving 
the time-dependent iV-particle Schrodinger equation. 
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Small Bose-Einstein condensates (BECs) of some 
1000 [1 or even 100 atoms [2] have been a topic of exper- 
imental research for several years. Recently, the investi- 
gation of many-particle wave- functions of BECs in phase 
space became experimentally feasible [3 j. This experi- 
mental technique will further investigations of beyond- 
mean-field (Gross-Pitaevskii) behaviour for small BECs. 

For a BEC initially loaded into one of the wells of a 
double-well potential, the many-particle oscillations of- 
ten seem to be damped compared to the mean-field be- 
haviour. Figure [I] shows such an apparent damping, 
which in fact is a collapse which will eventually be fol- 
lowed by at least partial revivals (cf. Refs. [H |5]), for 
N = 100 particles. This apparent damping coincides 
with an increase of the fluctuations of the number of par- 
ticles in each well [Fig. l](b)]. 

In order to numerically calculate the many-particle dy- 
namics, the Hamiltonian in the two- mode approxima- 
tion [6] is used, 



Hn = -J 



where are the boson creation and annihilation oper- 
ators on site j, hj = c^-c- are the number operators, J is 
the hopping matrix element and U the on-site interaction 
energy. 

The experimentally measurable [7] population imbal- 
ance is useful to quantify the oscillations depicted in 
Fig.0 
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FIG. 1. (Colour online) a) Population imbalance (z)/2 
[Eq. ([2])] as a function of dimensionless time r ^ for small 
BEC in a double-well potential. Initially, N — 100 atoms 
are in well 1, the quantum dynamics is given by the Hamilto- 
nian Q (NU/J = 0.4). b) Variance of the population imbal- 
ance Q . Experimentally, it will be difficult to distinguish the 
apparent damping and the increased fluctuations (which are 
both triggered by a collapse and revival phenomenon) from 
true damping introduced, e.g., by decoherence. 



The variance of the population imbalance can be quanti- 
fied by using the experimentally measurable [7 quantity 
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with < F z < N. For pure states, Eq. Q coincides with 
a quantum Fisher information [8 . Like the spin-squeezed 
states investigated in Ref. [7 (and references therein), 
quantum superpositions with large fluctuations are also 
relevant to improve interferometric measurements be- 
yond single-particle limits. A prominent example of a 
quantum superposition relevant for interferometry are 
the NOON-states [9] 
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(5) 



i.e., quantum superpositions of all particles either being 
in well one or in well two; |tii, 712) refers to the Fock state 
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with rt\ particles in well 1 and ri2 particles in well 2. 
Suggestions how such states can be obtained for ultra- 
cold atoms can be found in Refs. [5j [TOHTB] and refer- 
ences therein. For pure states, F z > 1 indicates that this 
quantum superposition is relevant for interferometry [8]. 
However, it remains to be shown that the increased fluc- 
tuations are really due to pure states rather than statis- 
tical mixtures. 

It might sound tempting to use the revivals investi- 
gated in Refs. [4] [5] to identify pure quantum states. 
However, while such revivals can be observed, e.g., for 
two-particle systems [17 , the situation for a BEC in 
a double well is more complicated. In principle, very 
good revivals of the initial wave-function should occur as 
long as the system is described by the Hamiltonian ([!]). 
While partial revivals can easily be observed, (nearly) 
perfect revivals might occur for times well beyond exper- 
imental time-scales - in particular if the experiment is 
performed under realistic conditions subject to decoher- 
ence effects [18]. It is thus not obvious how such an ap- 
parent damping might be distinguished experimentally 
from decoherence effects which would lead to statistical 
mixtures with (now truly) damped oscillation similar to 
Fig. [I] The focus of this paper thus lies on an experimen- 
tally realisable "echo" technique to distinguish statistical 
mixtures from quantum superpositions by using periodic 
shaking. 

Periodic shaking [19] is currently being established ex- 
perimentally to control tunnelling of BECs [20 -25 . For 
the model 0, periodic shaking can be included via 
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where K is the strength of shaking and uj its (angular) 
frequency. For large shaking frequencies [26] and not-too- 
large interactions, the time-dependent Hamiltonian (|6| 
can be replaced by a time-independent effective Hamil- 
tonian: 

U 2 

#eff = - Jeff (c\c 2 + 4&L ) + — ^ ™j ~ *) ( 7 ) 



with 
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where Jq is the Bessel-function depicted in Fig. [2] (b). 
Such effective Hamiltonians have been successfully tested 
experimentally in optical lattices, see, e.g., Refs. [20] [27]; 
negative J e ff have been experimentally investigated in 
Refs. [22j[^]. There are, however, also examples [281129] 
for which two or more Bessel function are needed to un- 
derstand the tunnelling dynamics. 

In the present situation, the effective description ^ 
offers the possibility to quasi-instantaneously switch the 
sign of both the kinetic energy (via shaking, cf. Fig. [5| 
and the interaction (via a Feshbach-resonance [30]). Con- 
trary to special cases where the wave-function [3TJ [32] 
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FIG. 2. (Colour online) a) Sketch of a double well which 
is shaken periodically to control tunnelling for ultra-cold 
atoms, b) For high shaking frequencies, the tunnelling rate 
is modified by the Jo-Bessel function [see Eq. Q]. The two 
squares indicate a pair of shaking amplitudes for which the 
Bessel function has equal modulus and opposite sign (x\ = 
1.691732695 and x 2 = 3.831705970 with \ J (x h2 )\ ^ 0.403 - 
it will be shown in Fig.[5](d) that it is not essential to precisely 
use these values). 



can be changed to obtain time-reversal, for periodically 
driven systems the Hamiltonian can be changed by quasi- 
instantaneously changing both the tunnelling term [by 
switching the shaking amplitude, e.g., between values 
shown in Fig. [5] (b)] and the sign of the interaction via a 
Feshbach-resonance [30] ; 



Hi. 



deal 



= / +if € ff(r = 0) : t < t 
-H eS (T = 0) : t>t 



(9) 



The corresponding unitary time-evolution is given by 



U(0,t) 




Hsh(t=0) \ . 
U J • 



r < r 
r > r 
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with perfect return to the initial state at r = 2tq. How- 
ever, the turning point ro has to be chosen with care; 
only by taking ro close to the maximum of the shaking 
can unwanted excitations be excluded (cf. Refs. [33H35] ). 
Recent related investigations of the influence of the ini- 
tial phase of the driving [replacing cos(a;t) in the Hamil- 
tonian ([6| by cos(cjt + (j>)\ can be found in Refs. [36H38] . 

In the following, the time-reversal is demonstrated by 
numerically solving the full, time-dependent Hamilto- 
nian |6| corresponding to the ideal time-reversal Hamil- 
tonian ([9| using the Shampine- Gordon routine [39] . 
Contrary to time-reversal schemes on the level of the 
Gross-Pitaevskii equation [40] [41] , here time-reversal is 
used to distinguish interesting quantum superpositions 
from statistical mixtures. Before implementing the time- 
reversal, Fig. [3] shows the wave-function for N = 100 
particles which were initially in one well. After several 
oscillations, the wave-function no longer is in a product 
state. Both the population imbalance and the phase can 
be measured experimentally [7j; in Fig. [3] the squared 
modulus of the scalar product with the atomic coherent 
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FIG. 3. (Colour online) Wave- function after apparent damp- 
ing dynamics versus typical product state a) Parameters as 
in Fig. [I] except for K = 1.691732695 and hw = 32 J. The 
wave-function at r = 307T is displayed as a function of pop- 
ulation imbalance z and phase 4> between the two wells [cf. 



Eq. (11)]. This quantum superposition can be characterised 
by F z ~ 38.8 and it could thus be used to improve interfer- 
ometric measurements, b) If all 100 atoms occupy the same 
single particle state (here: z = 0, cj) = 7r/2), the wave- function 
is much narrower and (as no product state satisfies F z > 1) 
would not be interesting for interferometry. 
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is plotted. The angle corresponds to a population im- 
balance of 
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Ideally, it should be possible to show that the wave- 
function of Fig. [3] (a) indeed is a quantum superposition 
by using the time-reversal of Eq. ([9| and than showing 
that 



(^end) = (z)(2r ) 



(13) 



is one: There is only one many-particle wave-function for 
which this is the case. Furthermore, the unitary evolution 
of solutions of the Schrodinger equation guarantees that 
for two different solutions |^i(2to)) = ?7(to, 2to)|'0i(to)) 
and \ip2 (2to)) = U(tq, 2ro)\ip2 (to)), the scalar product 
would be the same at r = To and at r = 2to (as WU = 

i). 

However, the Hamiltonian ([9| is a high-frequency ap- 
proximation and it has thus to be shown that this works 
for realistic driving frequencies (cf. Fig. [4|. Further- 
more, although there is only one wave- function at r = To 
which exactly leads to the value (z) en( ± = 1 at r = 2to, 
other (less interesting) wave-functions might lead to val- 
ues close to (2;) en d = 1 (cf. Fig. [5]). 

Figure [4] shows that the time-reversal dynamics is in- 
deed feasible. On time-scales for which there is not even a 




-0.5 tl_ 



1000 



2000 




0.45 - 



FIG. 4. (Colour online) Time-reversal of the quantum dy- 
namics of a small Bose-Einstein condensate in a periodically 
shaken double well [Eq. (|6|]. a) Population imbalance |2| 
as a function of time r = tJ/h for the same parameters as 
in Fig. [3] (b) b) Red/dark solid line: all other parameters 
as in panel (a) except for r > to: Kq — 3.831705970 J and 
U = —0AJ/N; light blue/grey dash-dotted line: as red/dark 
line but hu = 12 J; in both cases the revival of the initial state 
is visible near r ~ 85. c) Population imbalance for the same 
situation as in panel (a) but for much longer time-scales, d) 
If the switching takes place continuously rather than instan- 
taneously [Eq. ( |14| ], the revival of the initial state can still be 
observed [same parameters as for the red / dark curve in panel 
(b)] (7 = corresponds to instantaneous switching). 



partial revival of the initial state characterised by (z) = 1, 
the proposed time-reversal dynamics leads to final values 
above (z en( ±)/2 = 0.45 (Fig. [5] shows that this is enough 
to show that the wave-function at r = To was indeed 
a quantum superposition). In order to show that the 
scheme does not rely on the switching to be truly instan- 
taneous at t = to [where to is linked to To via Eq. ([3|], 
the amplitude in Fig [4] (d) was switched according to 



1 + tanh 



(14) 



the switching between the two interaction values was cho- 
sen analogously (for instantaneous switching, the switch- 
ing time can also slightly deviate from the ideal switching 
time, it just has to be close to the shaking maximum). 

Figure [5] primarily investigates how close the value of 
(z en d) has to be to one in order to show that the inter- 
mediate state shown in Fig. [3] (a) really is a quantum 
superposition. Figure [5] investigates the time-dynamics 
of product states ( pTj ) under the time-evolution which for 
the quantum superposition of Fig. [3] (a) would lead to a 
revival of the initial state [Fig. [4] (b)] . Figure [5] (a) shows 
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on short time-scales, as can be seen in Fig. [5] (b) which 
uses a (z m in) which is analogously defined to Eq. (15) to 
calculate 
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FIG. 5. (Colour online) a) Two-dimensional projection of 
(zmax)/2 + 0.5 as a function of both the initial phase and the 
initial population imbalance [cf. Eq. (11); zo = cos(#o)] For 
the numerics, the Hamiltonian leading to the red/dark curve 
in Fig. El (b) for r > To was used, b) Two-dimensional projec- 
tion of Az [Eq. ( |16| ] shows that (z) does not change dramati- 
cally on short time-scales, c) F z (tq) Q as a two-dimensional 
projection, d) If the time-reversal scheme of Fig.[4](b) for the 
red/dark curve is repeated with not ideal driving amplitudes, 
(^end)/2 + 0.5 (shown as a two-dimensional projection as a 
function of both driving amplitudes normalised by their ideal 
values - cf. Fig. [2} still lies well above the values shown in 
panel (a). 



In addition to not approaching (z) = 1, many product 
states lead to very large fluctuations [Fig. [5] (c)]; these 
fluctuations are particularly large if one compares them 
with the tiny values of F z (2tq) c± 0.4 for the red/dark 
curve in Fig. [4] (b). This offers an additional route to 
distinguish quantum superpositions as in Fig. [3] (a) from 
statistical mixtures by carefully investigating how the 
product states dTT] ) with large contributions to Fig. [3] (a) 
behave. Figure p](d) shows that the time-reversal scheme 
is feasible even if the driving amplitudes only approxi- 
mately meet the ideal values [Fig. [5] (d)]. 

To conclude, time-reversal via quasi-instantaneously 
changing the sign of the effective Hamiltonian is exper- 
imentally feasible for ultra-cold atoms in a periodically 
shaken double well. The change of the sign of the Hamil- 
tonian is achieved by changing both the driving ampli- 
tude and the sign of the interaction; a particularly useful 
initial state is the state with all particles in one well. The 
numeric investigations show that the revival of the initial 
state can be used to distinguish damping introduced via 
decoherence from the apparent damping related to a col- 
lapse phenomenon. Even if the revival of the initial state 
is not perfect, the scheme clearly distinguishes product 
states from quantum superpositions with potential inter- 
ferometric applications. 



that 

(z max ) = max{(z)(r)}| 099ro < r < ro (15) 

lies well below the values achieved in time-reversal 
[Fig.[4](b)]. Furthermore, it does not change dramatically 
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